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Limit  Cycles  of  Planar  Quadratic  Differential  Equations 


D.  E.  Koditschek  and  K.  S.  Narendra 


Introduction: 

Since  Hilbert  posed  the  problem  of  systematically  counting  and  locating  the 
limit  cycles  of  polynomial  systems  on  the  plane  in  1900,  much  effort  has  been 
expended  in  its  Investigation.  A  large  body  of  literature  —  chiefly  by  Chinese 
and  Soviet  authors  —  has  addressed  this  question  in  the  context  of  differential 

his  paper 


equations  whose  field  is  specified  by  quadratic  polynomials, 
w ^ consider  the  class  of  quadratic  differential  equations  which  admit  a 


unique  equilibrium  state,  and  establish  sufficient  conditions  for  the  existence 
and  uniqueness  of  limit  cycles.  The  work  is  based  upon  insights  and  techniques 
developed  in  an  earlier  analysis  of  such  systems  p^motlvated  by  questions  from 


mathematical  control  theory. 


V 


Until  the  fifties,  work  on  quadratic  systems  chiefly  concerned  the  existence 
of  a  center.  In  1952,  Bautin  [2]  showed  that  a  given  equilibrium  state  can 
support  as  many  as  but  no  more  than  three  limit  cycles  under  a  quadratic  field. 

Three  years  later,  a  paper  by  Petrovskil  and  Landis  [3]  purported  to  show  that  a 
quadratic  system  could  support  no  more  than  three  cycles  on  the  entire  plane. 
Although  this  result  was  called  into  question  by  several  researchers  (and  the 
authors  later  acknowledged  an  error  in  the  proof  [4])  it  apparently  inspired  a 
number  of  attempts  over  the  next  decade  to  complete  the  Hilbert  program  for  quad¬ 
ratic  differential  equations  [5,6,7].  Since  Coppel's  useful  survey  [8],  to  our 
knowledge  there  have  been  only  two  contributions  to  this  problem.  Perko  [9]  has 
demonstrated  the  existence  of  limit  cycles  for  a  certain  class  of  quadratic  sys¬ 
tems,  and,  recently,  Shi  Songling  [10]  has  presented  a  quadratic  system  which  appar¬ 
ently  has  four  limit  cycles,  settling  the  question  of  the  validity  of  [3]  negatively 
Thus,  the  Sixteenth  Hilbert  Problem  remains  unsolved  even  for  quadratic  systems. 


Al 


By  a  quadratic  system"  Is  meant  the  differential  equation 


Ax  + 


T 

x  Gx 


2x2  T 

where  A,G,H  e  R  (and  x  Gx  denotes  the  scalar  product  of  the  vectors  x  and 
2 

Gx  e  R  ) .  He  adopt  the  convention 


B(x)  ^ 


T  "1 
xaGx 

T 

_  x  Hx 


Our  results  may  be  sumnarized  briefly  as  follows.  In  section  2  it  is  shown  (Lemma  1) 
that  any  quadratic  system  with  a  unique  equilibrium  state  may  be  written  in  the 


x  ■  Ax  +  c  xDx  (2) 

2  2x2 

where  c  e  R  ,  D  e  R  ,  and  (Lemma  2)  that  the  global  boundedness  of  all  solutions 
depends  upon  the  spectrum  of  the  pencil  A  +  VD,  when  p  takes  values  in  the  range 
of  a  specified  nonlinear  functional.  In  section  3,  Theorem  2  establishes  the 
existence  of  limit  cycles  based  upon  the  spectrum  of  A  and  A  +  pD  as  an  easy 
extension  of  the  previous  results.  Theorem  3,  in  section  3,  uses  a  functional 
specifying  the  eigenvalues  of  A  +  pD  to  demonstrate  that  such  systems  have  only 
one  limit  cycle. 

2.  Preliminary  Discussion  of  Techniques: 

As  mentioned  earlier ,  the  techniques  of  [1]  are  used  extensively  in  this 

paper,  and  are  introduced  here  for  ease  of  exposition.  We  will  use  the  follow- 

2  T 

lng  notation  throughout  the  paper:  if  x,  y  e  R  ,  then  x  y  denotes  the  scalar 
product  of  x  and  y;  |x,y|  denotes  the  determinant  of  the  array  formed  by  the 

fo  -ll 

coordinates  of  x  and  y.  The  skew  symmetric  matrix  J  ■*  I  I  represents  a  90° 

L1  qJ 

by  x  *  Jx. 


rotation  and  the  orthogonal  complement  of  x  is  denoted 


We  note  that 


Vi) 


£  '  _  o' 

2  k  r  :  4J 

•  wo 
•.  *-•  **  m  m 


i  hr  s 
°  ^ 
>  ►  g 

AlSi  < 
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|x,y|  -  yTxA  -  yTJx. 


A  line  through  the  origin  containing  x  is  denoted  <x>  ■  {y:y  •  ax,  a  e  R}. 


2x2  T 

If  A  e  R  ,  the  symmetric  part  of  A  is  denoted  Ag  “  1/2[A  +  A  ],  If  A  is  singular 


T  2 

it  can  be  written  as  the  outer  product  A  ■  ab  where  a,  b  e  R  .  The  zeros  of  the 


T  T  2 

quadratic  form  x  Ax  ■  x  A  x  then  lie  on  <aA>u<bx>.  Conversely,  if  for  some  a  e  R  , 


T  T  2 

x  Ax  is  zero  on  <*x>  then  A  “  [ab  ]  ,  for  some  b  e  R  . 


By  limiting  attention  to  quadratic  systems  with  a  single  equilibrium  state. 


the  structure  of  (1)  is  greatly  simplified.  This  is  a  consequence  of  the  fact 

2 


that  A  and  B  are  homogeneous  polynomials  in  x,  hence  map  subspaces  of  R  into 
subspaces. 


Lemma_l:  If  A  is  bijective  then  system  (1)  has  a  unique  equilibrium  state  only 


2  2x2  T 

if  there  exists  a  c  c  R  and  a  D  e  R  such  that  B(x)  *  c  xDx. 


Proof;  If  for  some  1  e  R  and  x^  e  R  ,  lAx^  “  B(Xq)  •  0  then  -  1/X  x^  is  an 
equilibrium  state  of  system  (1).  Hence,  we  require  that  I  Ax,  B(x)|  *  0  imply 


B(x)  ■  0  for  any  x  ^  0.  Since  |Ax,B(x)|  •  x^  q(x2/xj) ,  where  q  is  a  cubic 


polynomial  in  x2/xp  there  exists  at  least  one  real  zero  of  q,  (say  vQ) .  Hence 

the  system  has  a  unique  equilibrium  state  only  if  B  ■  0  on  <[  v  ]> .  This,  in  turn, 

®  T 

implies  that  both  quadratic  forms  in  B  share  a  common  zero  line,  or  G  «[cd  ]  , 

j  T  is 

T  i  2  al 

H  -  [cd»  1  where  c  -  [  x  ]  and  d, ,<L  e  R  .  Defining  D  -  [  _  ] ,  we  have  the 

*■  S  u  1  Z  j  T 


desired  result. 


0  1 


□ 

Since  it  has  been  shown  [8]  that  no  limit  cycles  of  (1)  may  occur  unless  the 
linear  part  of  the  field  has  complex  conjugate  eigenvalues.  Lemma  1  achieves  an 
immediate  simplification  of  the  problem.  Any  quadratic  system  with  a  unique  equi¬ 
librium  state  capable  of  supporting  a  limit  cycle  must  admit  the  expression  given 
in  (2). 

Representing  a  quadratic  system  in  the  form  (2)  affords  analysis  by  well  known 
methods  of  linear  algebra.  The  following  lemma  relates  the  properties  of  the 
pencil  A  +  MD  to  the  Interaction  of  the  quadratic  and  linear  part  of  the  field  in 
(2),  providing  the  basis  for  much  of  our  global  analysis. 


-4- 


2:  For  any  matrices  A,D  e  R2*2  and  any  x  e  R 1  %  if  y(x)  -  - 


Dx,x 


Is 


defined,  then  x  la  an  eigenvector  of  the  pencil  A  +  y(x)D  with  corresponding 


Proof;  Define  <*(x)  ^  |Ax,x|  and  6(x)  A  |dx,x| .  Since  |[6A-aD]x,x|  ■  6|Ax,x| 


2 

-  a|Dx,x|  -  0  for  all  x  e  R  »  It  follows  £hat  [6^-aD]x  *  n(x)x  for  some  real  valued 


T  T 

function  n«  But  x  "xn  ■  x  [ 5 A-aD] x  ■ 


x  Ax  x  Dx 


T  T 

x  JAx  s  JDx 


[x,J  x][Ax,Dx]|  “  -Xx|Ax,Dx| 
□ 


Hence  n(x)  -  -| Ax,Dx|  and  the  result  follows. 

Corollary  2.1:  Let  A  have  complex-conjugate  eigenvalues.  If  the  pencil  A  +  uD 
Is  bounded  and  non-singular  then  system  (2)  has  a  unique  equilibrium  state  at  the 
origin. 

Proof:  If  A  has  no  real  eigenvectors  then  p(x)  cannot  be  bounded  unless  |Dx,x| 

Is  sign  definite  -  i.e.  D  has  no  real  eigenvectors.  According  to  Leona  2 
A  +  |iD  Is  nonsingular  If  X(x)  +  0;  since  | Dx,x|  Is  sign  definite,  this  is  equivalent 
to  the  condition  |Ax,Dx|  is  sign  definite.  The  latter  assures  that  |Ax,B(x)|  “  0 


only  on  <c^>  where  B  s  0.  Thus  the  origin  Is  the  unique  equilibrium  state.  □ 


Corollary  2.2:  Let  A  have  complex-conjugate  eigenvalues.  If  the  pencil  A  +  wD 
Is  bounded  and  has  negative  real  eigenvalues  over  the  range  of  y(x)  then  system  (2) 
has  no  unbounded  solutions. 

Proof:  As  In  the  previous  corollary  the  fact  that  y  Is  bounded  immediately  implies 
that  D  has  complex-conjugate  eigenvalues.  If,  additionally,  the  pencil  has  a 
spectrum  on  the  left  helf  of  the  real  line  In  C,  then  the  sign  of  |Ax,Dx|  must 
agree  with  the  sign  definite  form  |Dx,x|.  We  will  now  Interpret  this  sign  agreement 


geometrically  to  show  that  the  spiral  curve  defined  by  a  single  loop  of  the  linear 
tD_ 


trajectory  e  y defines  a  poeltlve-lnverlant  region  In  the  phase  plane,  arbitrarily 
far  from  the  origin. 


Choose  a  point,  y,  on  <c>  whose  sign  is  opposite  to  the  sign  of  the  real 
part  of  the  eigenvalues  of  D,  say  on  the  positive  ray.  Let  A  ■  {etDy|t  e  10, t  ]; 
e*  Dy  ■  yy;  0  <  Y  <  1}  be  a  complete  spiral  loop  and  let  A  *  {?y|c  e  [y,l]}  join  its 
end-points  as  depicted  in  Figure  1. 

The  normal  to  the  curve  at  any 
point  x  e  A  lies  in  <JDx>  and  since 
x  JDx  ■  |Dx,x|,JDx  is  either  interior 
directed  or  exterior  directed,  de¬ 
pending  upon  whether  |Dx,x|  Is 
negative  or  positive,  respectively. 

With  no  loss  of  generality,  we  assume 
|Dx,x|  <  0,  hence  JDx  is  the  interior 
directed  normal  to  A  at  x.  Similarly, 

Jy  is  the  interior  directed  normal  to 

A  for  any  y  e  A.  It  now  suffices  to 

Y  Figure  1 

show  that  f  (x)JDx  >  0  for  x  e  A, 

T 

and  f  (y)Jy  >  0  for  y  e  A. 

Expanding  the  first  inequality,  we  have 

W  W  fp  ip  fp  W  tp 

f  JDx  ■  x  [A  +  c  xD  ]JDx  ■  x  A  JDx 
-  |Dx,Ax|  -  -|Ax,Dx|  >  0 

2 

for  all  x  e  R  .  Expanding  the  second  inequality,  we  have 

fTJy  »  -yTJf  ■  -yTJAy  -  cTy  yTJDy 
T 

-  -|Ay,y|  -  c  y|Dy,y| 

T 

hence,  because  c  y  >  0  for  y  e  A,  and  |Dy,y|  <  0,  the  desired  Inequality  holds  when 
the  second  term  dominates  the  first  term  far  enough  away  from  the  origin. 


—6— 


Corollary  2.3:  Let  A  have  complex  conjugate  eigenvalues.  If  the  pencil  A  +  yD 
Is  bounded  and  has  positive  real  eigenvalues  over  the  range  of  ii(x) ,  then  the 
system  (2)  has  unbounded  solutions  for  every  initial  condition  outside  a  compact 
neighborhood  of  the  origin. 

Proof:  Let  y  be  a  point  on  <c>  whose  sign  is  the  same  as  the  sign  of  the  real 

part  of  the  eigenvalues  of  D,  say  on  the  negative  ray.  Let  A  and  A  be  as  in 

the  proof  of  corollary  2.2,  depicted  in  Figure  1.  Assume  again  with  no  loss  of 

generality  that  JDx  is  the  interior  directed  normal  to  A  at  x  and  Jy  the  interior 

T  T 

normal  to  A  for  y  e  A.  We  need  to  show  that  f  Dx  <  0  for  x  e  A  and  f  Jy  <  0  for 

y  e  A.  Since  |Ax,Dx|  has  opposite  sign  to  |Dx,x|  under  the  assumption  that  the 

pencil  has  positive  real  eigenvalues,  the  first  inequality  follows  for  every  spiral 

loop  A.  The  second  inequality  holds  on  A  outside  of  the  last  loop  for  which 

T 

7  JAy 

“T - 

y  JDy 

The  reader  should  observe  that  the  conditions  of  these  corollaries  may  be 

tested  using  standard  computations  involving  the  elements  of  A  and  D  as  established 

2x2 

in  the  proofs.  To  summarise,  the  condition  that  a  matrix,  H  e  It  »  have  complex 

conjugate  eigenvalues  is  equivalent  to  the  condition  that  |Mx,x|  never  vanish  on 
2 

R  ,  or,  equivalently,  that  lJM]g  be  a  positive  or  negative  definite  matrix.  If 
A  has  complex  conjugate  eigenvalues,  then  w(x)  is  bounded  if  and  only  if  |Dx,x| 
has  no  zeros,  or,  equivalently,  D  has  complex  conjugate  eigenvalues  as  well. 

Under  these  conditions  the  sign  definiteness  of  A(x)  depends  entirely  upon  the 
sign  definiteness  of  |Ax,Dx|  -  |d|*|d  *Ax,x| ,  hence,  upon  whether  D  *A  has  complex 

2x2  “1  I  I  j 

conjugate  eigenvalues,  or,  using  the  identity  in  tC  ,  D  ■  J  D  J,  whether  [D  JA]g 
is  a  positive  or. negative  definite  matrix. 

It  is  worth  noting,  in  passing,  that  the  conditions  for  the  boundedness  of 
system  (2)  resulting  from  Lemma  2  and  its  corollaries  may  be  extended  to  cover 


T  I 

|yc  y|  is  less  than  the  constant 


any  planar  quadratic  8yatem  (1).  Earlier  work  [1], [11], [12]  has  established 
that  homogeneous  quadratic  systems  which  may  not  be  written  In  the  form  c  xDx 
must  be  unstable.  It  may  be  shown  [1],[13],  In  consequence,  that  system  (1)  must 
have  unbounded  solutions  If  It  cannot  be  written  in  the  form  (2) .  Imposing  the 
added  condition  that  the  linear  part  of  the  field  not  be  unstable,  and  adjusting 
for  special  cases  permits  the  following  characterization  of  any  globally  asymp¬ 
totically  stable  quadratic  differential  equation. 

Theorem  1  CKodltschek  and  Narendra  fill:  System  (1)  is  globally  asymptotically 
stable  if  and  only  If 


(1)  The  eigenvalues  of  A  have  non-positive  real  part; 

2  2x2  T 

(il)  There  exist  a  c  e  R  and  D  e  R  such  that  B(x)  *  c  xDx; 


(ill)  The  pencil  A  +  y(x)D,  where  y(x)  *  ' VJ** 


,  has  non-positive 


eigenvalues  for  all  x  £  R  ;  is  unbounded,  if  ever,  only  on  an  eigenvector  of  A 
in  the  null  space  of  D  with  | A |  +  0;  is  singular  on  at  most  a  unique  line  which 
coincides  with  <c^>  iff  | A f  i*  0. 

In  fact,  according  to  [13  ],  conditions  (ii)  and  (ill)  of  this  theorem  (if  we  relax 
the  stipulation  that  |a|  +  0  when  the  pencil  is  unbounded)  are  necessary  and  sufficient 
for  the  boundedness  of  solutions  to  any  quadratic  system  (1)  as  well. 

However,  in  the  sequel,  we  will  confine  our  attention  to  quadratic  systems 
of  the  form  (2),  and,  specifically,  those  guaranteed  by  Corollary  2.1  to  have  a 
single  equilibrium  state. 

Existence  and  Uniqueness  of  Limit  Cycles: 


He  now  apply  the  techniques  developed  in  the  previous  section  to  the  subject 
proper  -  an  account  of  the  limit  cycles  of  system  (2) .  While  those  results  will 
be  seen  to  establish  the  existence  conditions  directly,  the  proof  of  uniqueness 
Involves  further  development. 

According  to  the  results  of  Lyapunov,  the  local  stability  behavior  of  system 
(2)  is  entirely  determined  by  the  linear  part  of  the  field  when  the  matrix  A  has 


eigenvalues  with  non-zero  real  part.  On  the  other  hand,  the  global  arguments  of 
Corollaries  2.2  and  2.3  depend  upon  the  spectrum  of  the  pencil  A  +  yD.  Using  the 
latter  result  to  guarantee  boundedness,  and  assuming  local  instability  by  the 
former  argument,  we  take  advantage  of  the  special  nature  of  limit  sets  of  planar 
dynamical  systems  established  by  the  Poincare-Bendlxson  Theorem  to  conclude  that 
a  limit  cycle  must  exist. 

Theorem  2:  System  (2)  has  a  limit  cycle  if 

(i)  the  matrix  A  has  complex  conjugate  eigenvalues  with  non-zero  real  parts 
and  (ii)  the  pencil  A  +  WD  is  bounded  with  non-zero  eigenvalues  opposite  in  sign 
to  the  real  part  of  the  eigenvalues  of  A  over  the  range  of  y(x) . 

Proof:  Assume  that  (i)  holds,  and  the  eigenvalues  of  A  have  positive  real  parts. 

Then  the  origin  is  totally  unstable,  hence  for  some  positive  definite  symmetric 
2  T 

matrix,  P,  It  -  {x|x  Px  <  y)  for  any  y  >  0  is  a  positive  invariant  set  of  system  (2) . 

If  A  +  u(x)D  is  bounded  with  non-zero  eigenvalues,  then  the  origin  is  the  sole 

critical  point  of  system  (2),  according  to  Corollary  2.1.  By  Corollary  2.2,  if  the 

eigenvalues  of  A  +  y(x)D  are  negative,  then  all  solutions  of  (2)  are  bounded:  in 

particular,  the  Jordan  Curve  AuA  bounds  a  positive-invariant  set,  J,  containing  the 

origin.  Thus  J  -  (x  Px  <  y)  is  a  compact  positive-invariant  set,  free  of  critical 

points.  In  consequence  of  the  Poincare-Bendlxson  Theorem,  the  positive  limit  set 

T 

of  a  trajectory  in  J  -  (x  Px  <  y)  must  be  a  limit  cycle  [14,  p.  232,  Thm.  9.3]. 

If  the  eigenvalues  of  A  have  negative  real  part  and  the  eigenvalues  of  A  +  y(x)D 
are  positive,  en  Identical  argument  concerning  negative  limit  sets  using  Corollary 
2.3  will  establish  the  existence  of  a  limit  cycle.  □ 

While  the  question  of  necessity  is  not  formally  addressed  in  this  paper, 
it  is  useful  to  remark  upon  the  existence  of  limit  cycles  when  the  conditions  of 
Theorem  2  ere  not  met.  Asstiming  (li).  Condition  (i)  is  certainly  necessary  for 
systems  of  the  form  (2)  to  support  a  limit  cycle:  Conpel  [8]  h/  i  shown  that  A 
must  have  complex  conjugate  eigenvalues;  when  A  has  pur''1  /  L^glnary  eigenvalues 


and  A  +  y(x)D  has  negative  eigenvalues  Theorem  1  guarantees  global  asymptotic 

stability,  while  a  similar  argument  establishes  that  all  non-zero  solutions  of  (2) 

grow  without  bound  when  A  +  v(x)D  has  positive  eigenvalues  in  this  case.  If  (i) 

2 

holds  and  A  +  p(x)D  has  a  zero  eigenvalue  for  some  x  c  R  then  system  (2)  has  at 
least  one  critical  point  distinct  from  the  origin.  If  (i)  holds  and  A  +  y(x)D 
is  bounded  with  non-zero  eigenvalues  whose  signs  agree  with  the  real  part  of  o(A) 
then  system  (2)  is  either  globally  asymptotically  stable  according  to  Theorem  1,  or 
can  be  shown  to  admit  only  unbounded  non-zero  solutions.  However,  if  (1)  holds, 
the  spectrum  of  A  +  p(x)D  is  sign  definite,  but  p(x)  is  not  bounded,  then  D  has 
real  eigenvalues  and  the  possibility  of  a  limit  cycle  remains.  As  will  be  seen 
below,  there  is  good  reason  to  suspect  that  system  (2)  cannot  support  a  limit 
cycle  unless  D  has  complex  conjugate  eigenvalues.  If  true,  this  would  imply  that 
the  conditions  of  Theorem  2  are  both  necessary  and  sufficient  for  a  quadratic 
system  (1)  with  a  single  critical  point  to  support  a  limit  cycle. 

We  now  proceed  to  show  that  the  limit  cycle  established  by  Theorem  2  is  indeed 
unique.  Along  the  way  we  will  restate  the  conditions  of  that  theorem  (Lemma  3, 
below)  and  provide  a  better  intuitive  sense  of  the  mechanism  underlying  the  Isolated 
periodic  solution.  This  is  achieved  by  a  transformation  to  polar  coordinates. 


Assuming  A  has  complex  conjugate  eigenvalues  we  may  always  find  a  coordinate 
system  (under  linear  transformation  of  the  state)  such  that  A  •  01  +  u>J  -  where 


I  - 


1  0 
0  1 


,  J  - 


0  -1 

1  0 


,  and  o,  a)  e  ft  -  and  c  ■ 


Then,  defining  the 


A  2  2  1/2  A 

polar  coordinate  transformation  p  *  [x^  +  Xj  ]  ,  9  •  arc  tan  ^ equation 


(2)  may  be  written  as 


p  »  1/p  f  (x)x  •  p[o  +  Pd(6)J 


0  -  -  »  ♦  pd(6) 


(3) 


I 


where  d  and  d  are  functions  of  8  only  and  are  defined  by 

T  T  T 

j  /  A\  A  A  x  Dx  j/  A\  A  ^  x  D  Jx 
d(0)  -  cos  0  — -  ;  d(0)  *  cos  0  — ^ -  . 

XX  XX 

Under  the  assumptions  of  Theorem  2,  d  is  sign  definite  for  0  e  [—tt/2, tt/2]  and  we 
assume,  with  no  loss  of  generality,  that  sgn  m  ■  sgn  d  >  0.  We  define 
n(8)  *  od(0)  -  u>d(0)  and  assert  the  following. 

Lemma  3:  The  following  conditions  are  equivalent  to  those  stated  in  Theorem  2, 
and  hence  are  sufficient  for  the  existence  of  a  limit  cycle  of  system  (3) : 


either  (i)  a  >  0  and  n  <  0 
or  (ii)  o  <  0  and  n  >  0 


for  0  c  [“ir/2,n/2] 


Proof:  Since  A  ■  al  +  a>J,  condition  (i)  of  Theorem  2  is  equivalent  to  one  of  the 
sign  conditions  on  a.  From  Lemma  2,  the  eigenvalues  of  A  +  p(x)D  are  given  by 


A(x)  *  {dx^xJ^  ”  ~  |Dx  x|  ^0lx,Dxl  +  -  -3-  .  Thus,  for  6  e  [-tr/2,ir/2] , 

the  sign  conditionson  d  and  n  are  equivalent  to  condition  (il)  of  Theorem  2.  □ 

As  reported  in  [8],  limit  cycles  of  quadratic  differential  equations  enclose 

convex  regions,  hence,  any  periodic  solution  of  (3)  must  have  an  angular  derivative, 

*  a  2  - 

0,  of  constant  sign:  no  limit  cycle  may  leave  the  region  C  *  {x  e  R  |ut  +  pd  >  0} . 

Consider  x(t;pQ>,  a  trajectory  in  C  originating  at  p^  -  a  point  on  the  negative  X£- 

axis.  For  some  >  t^  >  0  we  must  have  x(t^;pQ)  ■  p  -  a  point  on  the  positive 

X2~axis;  and  xCtjjPg)  "  P£  “  a  point  on  the  negative  x^-axis  -  as  depicted  in 

Figure  2.  Denote  the  resulting  curve  in  the  right  half  plane  over  the  interval 

[0, txJ  as  rx,  and  the  left  half  plane  curve,  over  the  interval  as  T2* 

Evidently,  Tj  may  be  expressed  as  x(s;pj)  where  s  e  [O.tj”*.].  Now  map  into 

as  follows:  since  0  is  sign-definite,  for  every  t  e  [0,t^]  there  exists  a 

unique  s  e  [0,  and  C  >  0  such  that 


x(s;p,)  -  -cx(t;pj. 
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For  convenience  we  obeli  denote  points  on  the  right  hand  curve,  r^»  by  r(t),  and 
on  the  left  hand  curve,  Tj,  by  £(a),  letting  p  ■  ||r(t)|  and  A  -  ||  t(s)||  -  cp* 

The  chief  advantage  of  this  sap  is  the  Induced  functional  dependence  of  s  on  t, 
hence  the  ability  to  write  a  differential  equation  for  p  and  A  using  the  sane 
angular  Interval.  From  (.4)  and  the  above,  we  have,  for  fixed  initial  conditions. 


o  +  pd 

id  +  pd 


a  -  Ad 

id  -  Ad 


0  e[-ir/2,ir/2] . 


(4) 


The  restatement  of  Theorem  2  in  Lemma  3  lends  added  insight  into  the  mechanism  by 
which  x(t;pQ>  grows  and  decays  on  r^ur^*  Considering  condition  (i)  of  Lemma  3, 
since  od  >  0  on  [~i»/2,ir/2],  the  fact  that  n  <  0  necessitates  d  >  0  on  that  Interval. 

e 

Hence,  from  (4),  while  p  must  Increase  on  T.,  A  becomes  negative  when  T,  enters 

^  1  A 

the  region  V  &  {x  e  R^|x  <  -o  }  in  the  left  half  plane.  Moreover,  C  has  a 

x  x 

boundary,  3C,  in  the  left  half  plane  and  n  <  0  implies  3C  c  P  -  i.e.  that  certain 
trajectories  contained  in  C  must  enter  Since  4~  A  -*■-<*>  as  fc(s)  3C,  the 

growth  of  a  trajectory  on  1*^  is  countered  with  increasing  effect  on  a  portion  of 
resulting  in  a  limit  cycle.  Notice  that  if  D  has  real  eigenvalues  then  3  is 
no  longer  sign  definite,  hence  d  may  not  be  sign  definite,  and  these  remarks  are  no 
longer  valid,  underscoring  the  importance  of  the  requirement  that  D  have  complex 
conjugate  eigenvalues. 

The  differential  equations  in  (4)  define  two  families  of  functions,  p(8;p^)  and 
A(0;Aq),  parametrized  by  initial  condition  on  the  negative  and  positive  X2~axes, 
respectively.  Observing  that  Aq  "  p(ir/2;p^)  and  that  the  vector  fields  in  (4)  are 
smooth  when  x  e  C,  we  may  explicitly  regard  p  and  A  as  functions  of  0  and  p^, 
continuously  differentiable  in  both  arguments.  Since  distinct  Integral  curves  of 
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autonomous  systems  defined  by  smooth  fields  remain  distinct  over  all  times,  we  have 

-JL-  p  >  0  and  4 —  X  >  0  for  all  8  e[-ir/2,ir/2] .  Hence,  the  function 
3p0  3p0 

A  X(ir/2»p0> 

*(P0>  "  p(-*72,p0)  ’ 

which  represents  the  ratio  of  the  magnitudes  of  the  end  points  of  the  curve 
(both  on  the  negative  x^axis)  Is  a  continuously  differentiable  function  of  pQ. 
Evidently,  is  the  integral  curve  of  a  limit  cycle  if  and  only  if  ♦  ■  1. 

The  proof  of  uniqueness  involves  a  demonstration  that  is  monotone  in  Pq  over 
an  interval  of  interest,  and  hence,  may  pass  through  1  at  most  once.  That  demonstra¬ 
tion  depends  upon  the  following  computation. 

Corollary  3.1:  Conditions  (i)  and  (ii)  of  Lemma  3  respectively  imply 

<»  <«$»«•*» 

<41>  >2'le  Cln  M  *  0/“1 

Pro of.  £">  -^-/4r-  *■£*-  1%-.  tvmw  v.  have 


38 


3P0  »P0 


383p0  3pq 


383p0  3p0 


£-0  <  -H  >  ■  1 1? ln  x  +  x 


— * - ?  1  and  TT  <  J  1  TZ  ln  p  +  p - rr  1  TT 

(w-Xd)2  3p0  3p0  30  30  (w+pd)  3p0 


Since  d  >  0,  we  have - — —  >  1,  and - <  1,  for  all  8  e[-*/2,*/2] 


ut  —  Xd 


u  +  pd 


Hence,  if  condition  (i)  of  Lemma  3  holds  then 


-  2 
(ctt-Xd) 


b> 


w  -  Xd 


Ii —  “  — —  ln  X  -  a/w 
38 


and 


p - V2  <  -J - "  |e  ln  p  "  o/“* 

(ttfpd)  (it  +  pd 


since  n  <  0  on  [-ir/2,ir/2]  and  substituting  from  C4). 


Since 


3p 


> 


0, 


this  implies 


^  Ci|,< «  m  X  -  o/B,  |1 


_L_  /£&>  <  (2 
3p0  '■ae 


36 


In  p  - 


o/w) 

3p0 


yielding  (i),  above.  The  identical  argument  holds  for  (li)  with  signs 
reversed,  since  n  >  0.  □ 

We  may  now  state  the  second  principal  result  of  this  paper. 

Theorem  3:  Under  the  conditions  of  Theorem  2,  system  (2)  has  only  one  limit 
cycle. 

Proof :  x(t;pg)  is  a  limit  cycle  of  (2)  if  and  only  if  KPq)  *  1  in  system  (3). 

According  to  Theorem  2,  L  ■  {p^  >  0|  4>(Pq)  "  1}  is  non-empty,  and  bounded  away  from 

the  origin,  hence  p*  ■  inf  L  exists  and  Pn  >0.  We  will  show  that  —■  «l>  is  sign  definite 

0  w  apn 

*  * 

for  all  Pq  >  Pq*  hence  x(t;p^)  is  the  only  limit  cycle  of  (2). 
j  la 

Note  that  — —  *  ■  “  [  — —  A  (w/2,PA)  -  ♦] .  We  will  show  below  that  under 
dp0  p0  3p0  0 
condition  (i)  of  Lemma  3, 


A  (v/2,Pq)  < 


* 


and  hence 


1]. 


Since  a/u  >  0  and  KPq)  ■  1»  this  is  clearly  negative  for  p^  >  p^.  Similarly, 

under  condition  (11)  of  Lemma  3  the  inequalities  are  reversed,  and  o/u  <  0 

d  * 

so  that  -r—  ♦  >  0  for  all  pA  >  pA. 

OPq  o  o 

3 

To  obtain  the  bound  on  — —  A  (it/ 2 

3P0 


,PQ)  we  recall  that  A(-»/2 ,Pq)  ■  P(tr/2,Pg. 


and  p(-»/2,Pq)  ■  Pq»  hence 
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»  (./2.P0)  »  C/2,P0)  -g-  p  C/2.P0)j-  X  (-./2,P0)^-<-,/2,, 

'  /  m  lB  [ir A  (e’co)  Ir 0  «’">o,ld9 


Applying  Corollary  3.1  to  case  CD  yields 


3  f  w/2  a 

^  3^  A  ("/2’p0>  <  j  2(jy  In  PA  -  o/o>)d0  -  In  ( 


A(ir/2,p  )  2 

7,  T7o  ?  -v  )  “  2*  o/w 


p(-ir/2,p0) 


In  ♦  -  2ir  o/«. 


^enCe  A  (ir/2»Pg)  <  ♦  e  claimed.  Case  (ii)  proceeds  identically  with 

the  signs  reversed.  q 

4.  Conclusions* 


This  paper  presents  sufficient  conditions  for  the  existence  of  limit  cycles 


of  quadratic  systems  with  a  unique  equilibrium  state.  The  conditions  guarantee 
that  the  limit  cycle  is  unique.  The  results  are  based  upon  insights  and  techniques 


developed  during  an  earlier  investigation  of  the  global  stability  properties  of  (1) 
ll],  facilitated  by  the  expression  of  that  system  in  the  form  (2).  They  strongly 


suggest  that  these  conditions  are  necessary  as  well,  hence,  that  no  quadratic 
system  with  a  unique  equilibrium  state  can  support  more  than  one  limit  cycle. 


That  result,  the  uniqueness  of  a  limit  cycle  around  any  equilibrium  state  of  (2), 
fhe  relation  of  limit  cycles  of  (2)  to  those  supported  by  general  quadratic  systems 
(1),  all  remain  to  be  rigorously  established. 
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